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Comment on "Violation of the Luttinger sum rule within the Hubbard model on a 
triangular lattice", by J. Kokalj and P. Prelovsek, Eur. Phys. J. B 63, 431 (2008) 
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Using the first-order series expansion of tfie function G(fc; y), in powers of /i = y — U /2, pertaining 
to the insulating ground state of a single-band Hubbard Hamiltonian at half-filling, Kokalj and 
Prelovsek (KP) have in a recent paper [Eur. Phys. J. B 63, 431 (2008)] reported breakdown of 
the Luttinger theorem for the specific case where the lattice on which the Hubbard Hamiltonian is 
defined is a two-dimensional triangular lattice, for which the ground state is not invariant under 
particle-hole transformation. Here G{k;fi) is the single-particle Green function G{k;e) evaluated 
at e — y, the zero-temperature limit of the chemical potential corresponding to half-filling, and U 
the on-site interaction energy. In this Comment we demonstrate that unless fl = (to be strictly 
distinguished from jl small but non- vanishing) , any finite- order series expansion for G{k; n) in powers 
of in general falsely signals breakdown of the Luttinger theorem. The violation of this theorem as 
asserted by KP is therefore an artifact of their first-order calculation. 
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Considering the non-magnetic uniform gromid state 
(GS) of Hamiltonian H, corresponding to N fermions, 
the associated single-particle Green function G°'{k;e) is 
independent of the spin index a. Consequently, in this 
Comment wc suppress a and denote the latter func- 
tion by G{k;e). For the mentioned iV-particle GS, 
the Luttinger theorem under consideration states that 
[1, 2, 3, 4, 5, 6, 7, 8] 



7V = 2^e(G(fc;/.)), 



(1) 



where ^ is the zero-temperature limit of the chemical 
potential satisfying the equation of state of the grand- 
canonical ensemble in which the mean number of parti- 
cles is equal to [7] . The summation with respect to k 
in Eq. (1) is over the entire wave- vector space available 
to the system under consideration. For instance, for H 
defined on a Bravais lattice, {Ri}, this space consists of 
the points constituting the corresponding first Brillouin 
zone, IBZ. 

In the thermodynamic limit, where the underlying k 
space consists of a continuum, one defines the Luttinger 
surface 5l corresponding to the insulating GS of H as 
the locus of the k points for which 



G(fc;/i) = 0. 



(2) 



We point out that Eq. (2) should be used with caution, 
as 5l has a more encompassing definition than implied 
by this equation (Sec. 2.4 in Ref. [7]). For the considera- 
tions in this Comment, we shall not strive generality and 
therefore in the following consider Eq. (2) as the defining 
equation for 5l. 

The Luttinger surface separates the k space into a re- 
gion whose points contribute to the Luttinger sum on the 
right-hand (RHS) of Eq. (1), i.e. the Luttinger sea, and a 
complementary region. For the two-dimensional system 
that we consider here, it turns out that the underlying 
5l consists of one piece, and so does the Luttinger sea. 



For the economy of notation, wc introduce the radial 
unit vector h centred at the origin of the fe space (here 
the IBZ) and directed along the radial vector k under 
consideration, whereby wc can write 



fc = fe n, 
where fc > 0. We thus define 

9" a" 



G'("'")(fc;/i) 



-G(fc;/i). 



(3) 



(4) 



With fl denoting the exact zcro-tcmpcraturc limit of 
the chemical potential corresponding to particles, for 
an insulating GS KP [9] employed the first-order approx- 
imation of the exact Taylor expansion 



^ — ^ n I 



ri=0 



where 



fi = fl - -U. 



(5) 



(6) 



The series in Eq. (5) being a power series, its validity is 
restricted to the region where |//| is less that the smallest 
of the two positive numbers £+ — [//2 and U/2 — in 
which e+ (e_) is the smallest (largest) energy e greater 
(less) than /i for which ^(fc; e) is non- vanishing; here 



A{k; e) = T- Im[G(fc; e ± iO+)] 



(7) 



is the single-particle spectral function. 

More explicitly, by employing the spectral represen- 
tation for G{k;ii) in terms of A{k;e), KP [9] expressed 
G{k; fj.) as 



G{k;fi) = G-{k;fi) + G+{k;fi), 



(8) 



2 



in which G^{k;fi) arise from the spectral contributions 
corresponding to £ ^ /x. Using the standard moments ex- 
pansion of the underlying integrals [10, 11], KP obtained 

that [9] 



ri=0 m=0 ^ ^ 

(9) 

where the explicit expressions for M^{k) are given in 
Eq. (6) and (7) of Ref. [9]. 
Making use of the identity 



(10) 



n^O rn— n— m—n 

the expressions in Eq. (9) can be cquivalently expressed 



G^ik■,^l) 



n=0 



m—n 



2 \ m-t-l / 



(11) 

from which the explicit expression for G^"'"^ (fc; U/2) can 
be immediately read off. We point out that since KP [9] 
considered the problem at hand in the regime \t/U\ ^ 1, 
in their explicit calculations they employed expressions 
for G(°^")(A;; U/2), n = 0, 1, that are correct to the lead- 
ing non-trivial order in t/U. This amounts to a specific 
truncation of the sum with respect to m on the RHS of 
Eq. (11). We shall briefly return to this aspect later. 

Denoting by G^(fc; /z) the function that one obtains on 
replacing the upper bound of the sum on the RHS of 
Eq. (5) by j/, one has 

— ^ n\ 

^G.(fc;/.)+ ^'"7'!^^;,^^^^ for A^O. (12) 

(i^ + 1)! 

Since G(fe; /i) is unknown, within the approximate 
framework adopted by KP [9] one calculates an approx- 
imate Luttinger surface, S^'^\ through solving the equa- 
tion 



G.ik;^l)^0, 



(13) 



where v = 1 for the explicit calculations reported by KP 
[9]. For a predetermined n, we denote the solution of 

Eq. (13) by fc['^\ and, with reference to the convention 

in Eq. (3), the corresponding k by ki'^'' . Similarly, for 
the solution of Eq. (2) in the direction of the same n, 
which we denote by fcL, and the corresponding k by kj^. 
Following Eqs. (12) and (13) one thus has 



Giki^'^-fi) 



G(0,.+l)(fcO 



(^ + 1)! 



for fi ^ 0. 



(14) 



Since kj^ as jl —>■ 0, by assuming that G(fe; ^i) is 

regular for fe in a finite neighbourhood of 5l, one can ex- 
pand G{k\^^ ; ^) in powers of k\^^ — ki^ for sufficiently small 
values of ji. Since by definition G{ki^\ /i) = 0, Eq. (2), one 
thus obtains that 



A J I 



(15) 



m— 1 



Truncating this series at m = 1, and assuming that 
G^-'-'^^ki^; fi) 7^ 0, from the resulting expression and the 
asymptotic result in Eq. (14) one deduces that 



G.(o,.+i)(fcM.^/2) 



G(i^")(fc,;;t/) (i^ + iy. 



for fl ^ 0. 



(16) 

For sufficiently small values of /i, one can replace the ki^ 
on the RHS of this expression by k['^^ , and vice versa. 

The simple asymptotic expression in Eq. (16) makes 
explicit that unless /i = 0, kj^"^"^ does not coincide with 
/cl for any finite value of i^. Evidently, it is in principle 
possible that k['^'' may be equal to /cl for some directions 
of n, however this equality cannot apply for all directions 
of n. Although in spite of /cl'^'' ^ fct, the Luttinger theo- 
rem, Eq. (1), may apply when the G(fc; fi) on the RHS of 
Eq. (1) is replaced by G,j{k; fi), with v < oo (see Fig. 47, 
p. 231, in Ref. [12]), this need not be the case in general. 

For instance, in the cases where both 5l and ^l"^' , with 
< oo, are isotropic or nearly isotropic, the Luttinger 
theorem in terms of Gi^{k; fi) unquestionably fails. We 
have thus unequivocally demonstrated that for v < oo, 
the strategy adopted by KP in Ref. [9] is not appropriate 
for verifying the validity of the Luttinger theorem. 

We shall now focus on the numerical results reported 
by KP [9]. As we have indicated earlier, these results 
correspond to u ~ 1. For the insulating GS of a single- 
band Hubbard Hamiltonian, in Ref. [9] (Eq. (16) herein) 
KP obtained that 



^(fc;/^) = ^(E^^^^-^^e^) -/^) (17) 
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where [9] 

Cs = {Ss ■ So) (18) 
is the GS spin-spin correlation function, and 



(19) 



in which (5 G {1, 2, ... } refers to nearest-neighbour, next- 
nearest-neighbour, . . . , so that ^ .^^ is a sum over the 
lattice vectors {Ris} which are nearest-neighbours, next- 
nearest- neighbours, . . . , to the central lattice vector 
for S = 1,2,..., respectively. With reference to the re- 
mark following Eq. (11) above, we note that the error 



3 



0{tij/U^) on the RHS of Eq. (17) corresponds to re- 
stricting the upper bound of the sum with respect to m 
on the RHS of Eq. (11) to 1. 

SpeciaUsing to the cases where the hopping integral ti^ 
is non- vanishing only for nearest-neighbour sites (assum- 
ing further that tij = t for Ri and Rj nearest neigh- 
bours) and suppressing the last term on the RHS of 
Eq. (17); KP [9] employed the following simple expression 
for t//t = 40 in examining the validity of the Luttinger 
theorem [13]: 

'(20) 

where [9] 



Ci = -0.182, 
and (identifying the lattice constant a with unity) 



(21) 



1 \/3 

ei(fc) = -2t(cos(fc^) -H2cos(-A:,0cos(— fcy)), (22) 

in which and ky are the Cartesian coordinates of k. 
With reference to Eq. (2), one observes that according 
to Eq. (20) for the system under consideration the iSl 

(i.e. the iSl^') corresponding to \U /t\ ^ 1 coincides with 
the Fermi surface of the non-interacting system, as- 
sociated with ei(fc) and corresponding to the chemical 
potential /j^o, where 



Mo 



4Ci 



(23) 



The value of fiQ corresponding to half-filling is obtained 
by solving the equation 



f d2fce(/io-£i(fc)) = i, 

IBZ JlBZ ^ 



(24) 



where J^ibz = Stt^/a/S is the area of the IBZ under 
consideration in the units where the lattice constant a 
is identified with imity; this IBZ is a regular hexagon 
[14, 15] of which the Cartesian coordinates of the ver- 
tices arc: 

27r(|0), 27r(i,±^), 27r(-i,±^), 27r(-|o). 

(25) 

Numerical calculation yields /ip ~ 0.852^, which following 
Eq. (23) leads to [9] 



/i = yui = 4CiMo ~ -0.620 1 



(26) 



as the value for /2 required for the Green function G(fc; fi) 
in Eq. (20) (i.e. G,y=i{k; /i)) to satisfy the Luttinger the- 
orem, Eq. (1). 

Through a combination of explicit numerical calcula- 
tions (see Fig. 1 in Ref. [9]), KP [9] deduced the following 
expression for p, specific to half-filled uniform GSs: 

6 8t^ 

fi (0.19±0.1)i+ for 0.025t < — < 0.05i. 



With U/t = 40; this expression yields 
fi^ (-0.02 ± 0.1) i. 



(28) 



U 



U 



(27) 



a range of values for /i which is far too removed from the 
value presented in Eq. (26), i.e. the value required in or- 
der for the Luttinger theorem to apply for large values of 
\U/t\, and sufficiently large systems (for which the sub- 
stitution of Y.k by ly/i2n)^]J^^^d'^k, where is the 
area of the system, is justified). Thus, KP [9] arrived at 
the conclusion that for the GS under consideration, the 
Luttinger theorem broke down. 

The fallacy in the reasoning by KP [9] becomes ap- 
parent by realising that the Green function in Eq. (20) 
concerns G^=i(fc;/x) (sec Eq. (12)), to be distinguished 
from the exact G{k;fi), for large values of \U/t\. With 
reference to Eq. (16), it is evident that unless the coeffi- 
cient function multiplying (1^^^ / {v-\-\)\ is vanishing (here 
v = 1), any non- vanishing value of fl implies deviation of 

the S'l'^^ calculated on the basis of G^{k\ /i) from the ex- 
act iSl- Since in the case at hand s'^^^ is nearly isotropic 
(see Fig. 2 in Ref. [9] — one expects the same to ap- 
ply for the exact S^), this deviation directly leads to an 
apparent; but false, violation of the Luttinger theorem. 

On general grounds, one can demonstrate that 
(^(0,1^+1) (;^(>'). fj/2) is for insulating GSs relatively large, 

leading to a considerable deviation of k'^^ from /cl even 
for relatively small but non- vanishing values of (see 
Eq. (16)). To appreciate this fact, one should realise that 
by the Kramers-Kronig relation for Re[G(fc; e)], the sud- 
den change of Im[G(fc; £±i0+)] for e at band edges (from 
identically vanishing for e inside the underlying gap, to a 
function whose magnitude is steeply increasing for e im- 
mediately past the band edges) implies very rapid change 
in Re[G(fc; e)] for e inside the gap region, leading to very 

large values of G(°''^+i)(fc!;''^t^/2) [16]. Below we specify 
the scale according to which the adjective 'large' is to be 
understood here. 

For a quantitative analysis, we note that for the coef- 
ficient of fi'^^^ /{v+ 1)! in Eq. (16), specific to v = 1, one 
expects that 

This result follows from the assumption that G{k;z) 
is a continuously diffcrentiable function of k and z in 
the neighbourhoods of 5l and U/2 respectively, whereby 
/cl^^ — fct must to leading order vary linearly with fl for 
small values of jl. On the basis of this observation, the 
relationship in Eq. (29) follows immediately on dividing 
both sides of Eq. (16) by ft (assuming that jl ^ 0). With 
reference to the result in Eq. (27) (or Eq. (28)), the es- 
timate in Eq. (29) implies that for the case at hand the 
left-hand side of this equation is to leading order propor- 
tional to 10/t. The scaling of the latter value with 1/t, 
to be contrasted with l/U, is most significant. 
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Following the observation that to leading order (/cl — 
k^'^)/ki^ is proportional to /2/t, from the expression in 
Eq. (28) one concludes that for the case at hand (/cl — 
k\^^)/ki^ must be of the order of 10%, which conforms 
with the value of 17% as reported by KP [9]. 

We thus conclude that the inference by KP [9], that 
the Luttinger theorem broke down for the particle-hole 
asymmetric insulating GSs of the single-band Hubbard 
Hamiltonian. is incorrect. The reported deviation of the 
calculated Luttinger surface from the Luttinger surface 
for which the Luttinger theorem would apply (for large 
values oi U/t), is accounted for by the truncation error 
associated with the use of the first-order expansion for 
G{k; fi) in powers of /i = fj, — U /2. The attempt by KP 
[9] at finding the exact solutions of Eq. (1) (the totality of 
which comprise the exact iSl) on the basis of the approx- 
imation G{k;fi) w Gi,=i{k] fi), which applies for a rela- 
tively small neighbourhood of /2 = 0, may be likened with 
that of finding the solutions of sin(x)/a; = at a; = ±7r 
on the basis of the approximation sin(a;)/a; « 1 — a;^/6, 
which applies for a small neighbourhood of a; = 0; the 
solutions of 1 — a:^/6 = are x = ±V6 « ±2.45, which 



are indeed considerably removed from ±7r (explicitly, by 
some 22%, which is comparable with the 17% reported 
by KP [9]). 

In view of the above observations, we wish to close this 
Comment by a critical message. As the detailed consid- 
erations in Ref. [7] have shown, the Luttinger theorem 
is valid under the conditions specified by Luttinger and 
Ward [1]. In Ref. [7] considerable amount of space was 
devoted to demonstrating that a vast body of the ex- 
tant publications, that supposedly demonstrate failure 
of the Luttinger theorem, are distinctly erroneous. This 
statement equally applies not only to Ref. [9] , considered 
here, but also to a most recent publication by KP [17] to 
which a separate Comment [18], by the present author 
and A. M. Tsvelik, is directed. To our best knowledge, 
to this date no claim of the supposed failure of the Lut- 
tinger theorem has withstood the test of time, so that 
it seems high time that henceforth researchers leave this 
theorem undisturbed, and resist the temptation of appar- 
ently interminably declaring it as invalid through mistak- 
ing the imperfections of their pertinent calculations with 
the failure of this perfectly valid theorem. □ 
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